FILE  COPY 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LINCOLN  LABORATORY 


MULTIPLE  BEAM  FORMING  NETWORKS 

W.  C.  CUMUINGS 
Group  61 


TECHNICAL  NOTE  1978-9 

18  APRIL  1978 


Approved  (or  public  releo^e;  distribution  unlimited. 


LEXINGTON 


MASSACHUSETTS 


ABSTRACT 


Lossless  beam  forming  networks  uclng  phase  shifters  and  quadrature 
couplers  arc  described.  A technique  for  using  che  Butler  network  with 
randomly  spaced  arrays  Is  discussed.  A synthesis  procedure  for  networks 
with  N-element  arrays  Is  described.  A method  of  using  these  networks 
with  planar  arrays  Is  discussed.  A least  mean  square  algorithm  for  che 
determination  of  a beam  forming  newtork  for  an  arbitrary  array  and 
arbitrary  beam  directions  is  developed.  Computed  pattern  data  for  two 
thinned  random  planar  arrays  Is  included. 
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INTRODUCTION 


Beam  formlrg  networks  are  generally  conoidered  to  be  circuits  which  couple 
two  or  more  antennas  to  one  or  more  input/output  terminals.  The  terminals  are 
considered  input  terminals  if  the  entennae  are  transmitting  antennas  and  output 
terminals  if  the  antennas  are  receiving  antennas.  A power  divider  network 
which  couplets  the  antenna  elements  of  a linear  array  to  a single  terminal  in 
such  a fashion  as  to  result  In  a linear  phase  taper  across  the  face  of  the 
array  is  one  fc  nn  of  c beam  forming  network.  With  this  network,  the  array 
produces  a bean  in  a direction  determined  by  the  array  element  spacing  and  the 
slope  of  the  phase  taper.  Similarly,  the  hybrid  which  is  comiected  to  two 
elements  to  form  the  monopulse  sum  anu  difference  patterns  is  a beam  forming 
netv’orkv 

Beam  forming  networks  can  be  linear  or  nonlinear,  reciprocal  or  nonreci- 
procal, passive  ut  active,  loeay  or  losslaas,  matched  or  not  matched.  They  can 
be  used  to  couple  many  antenna  elamenta  to  a single  port  or  c few  antenna 
elements  to  many  ports.  This  technical  note  is  concerned  with  a single  type  of 
beam  forming  network.  The  network  in  liuaar,  reciprocal,  passive,  loasless, 
matched  and  connects  N entanna  elements  to  N ports  to  fora  N antsnna  besms. 

The  network  discussed  utilizes  two  circuit  elements,  the  3-dB  quadrature 
coupler  and  the  fixed  phase  shift. 

The  discussion  will  coniine  itself  to  the  use  of  the  Individui-J.  antenuas 
as  slemento  of  an  array.  Initially,  soma  restrictions  will  be  placed  on  the 
antenna  array  configuration.  These  restrict Iona  will  ba  removed  as  the  dis- 
cussion prograases.  First,  the  Butler'*'  Matrix  is  briefly  dlacussad.  Next  the 
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Nolan  Matrix  is  described.  Methods  for  simpllfylnR  the  Nolan  Matrix  arc 
discussed,  and  the  Butler  Matrix  is  shown  tc  be  a special  case  of  the  Nolan 
Matrix.  Finally,  an  extension  of  ths  Nolan  Matrix  to  completely  arbitrary 
array  configurations  is  described.  Computed  radiation  pattern  data  arc  pre- 
sented for  some  arbitrary  arrays. 


II.  TWO  SIMPLE  BEAM  FORMING  NETWORKS 


The  simplest  form  of  the  N temlnal  by  N antenna  beam  forming  network  Is 
shown  In  Figure  1.  This  simply  connects  each  antenna  element  to  each  Input/ 
output  port,  and  each  beam  formed  is  the  beam  of  a single  element.  This  type 
of  network  may  appear  trivial  but  In  many  situations  it  can  be  the  best  solu- 
tion to  a design  problem.  The  beams  formed  are  completely  Independent,  the 
array  configuration  Is  arbitrary  and  not  relevant  to  the  beams,  and  failure  of 
a single  antenna  has  an  effect  on  only  one  beam. 

Many  design  prublens,  however,  require  that  the  be^ms  be  narrow  and  that 
the  aperture  size  of  the  antenna  array  be  small;  two  requirements  that  are  In 
conflict  with  one  another.  For  this  situation.  It  Is  desirable  to  form  each 
beam  using  as  much  of  the  available  aperture  space  as  is  possible.  Each  I/O 
port  should  couple  to  mors  than  one  of  the  sritenna  elements  In  the  array. 

Figure  2 Illustrates  a method  for  doing  this  with  a two-element  array. 

The  coupler  convention  utilized  is  indicated  to  the  right  of  the  figure.  Note 
that  un  input  to  port  I results  in  antenna  2 lagging  antenna  1 by  90*  in  phase 
causing  Che  beam  formed  to  squint  to  the  right  in  tho  figtvre  while  an  Input  to 
port  2 causes  the  beam  to  squint  to  Che  left.  The  amcunt  of  beam  squint  and 
tne  width  of  the  beame  formed  ie  piimerily  dependent  on  the  distance  between 
the  antennas.  Hie  beam  forming  natwork  doas  not  raqulre  any  specific  antenna 
element  spacing;  two  beairz  will  be  formed  regardless  of  the  spacing.  If  the 
spacing  Is  large  anough,  the  radiation  patterns  for  each  of  the  beeiss  will 
include  grating  Icbae.  If  the  antenna  requirement  is  such  that  It  is  desired 
to  point  the  two  beams  formed  nt  a narrow  field  of  the  view  and  the  pattern 
structure  outside  this  field  of  viow  is  of  no  consequence,  then  the  spacing  can 
be  increased  to  the  point  where  either  the  grating  Iches  bej^ln  to  illumlnar.e 
the  field  of  view  or  the  width  of  the  two  formed  beams  reaches  a desired 
minimum.  Such  a raqulisment  might  apply  for  a space  satellite  antenna  illu- 
miuatlng  the  disc  of  the  earth,  for  axample. 
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ANre.'4KA  2LCMCMTS 


■EAM  FONMINO  NETWORK 

Fig.  1.  N X N beaaforalng  nctvork-alBplwRt  fora. 


Fig.  2.  TWo  cleaentt  two  bona  array. 
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G«n«rally,  for  boat  parforvanca  of  tha  catwork  ahovn  in  Flgura  2,  tha 
antauna  olaaont  radiation  paf.tama  ahould  ba  aubatantially  tha  aaaa,  and 
tha  alaaanta  should  ba  pointad  in  tha  aaaa  dlractlon.  Obviously,  tha  t%ro- 
alanont  array  la  linaar.  Thasa  rastrlctions  ara  iapoaad  on  tha  arrays  which 
will  be  discusaod  in  tha  naxt  saction. 
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III.  THE  BUTLER  BEAM  FOBMIMG  NETWORK 


Butlar  Bean  Foralng  Networks  have  been  described  extensively  Jn  the 
13  4 

literature  * * . This  discussion  will  suaBsrlze  sons  of  thess  descriptions 
with  the  object  of  extending  the  bean  fomlng  natwork  concept  to  a class  of 
irregularly  apscad  llneor  arrays. 

Begin  with  the  array  restrictions  Indicated  In  the  previous  section  and 
further  that  tha  elanants  In  the  linear  array  are  equally  spaced  and  that  their 
number  Is  equal  to  an  intagral  power  of  two.  Since  there  are  an  even  number  of 
elements,  hence  beams.  It  la  reasonable  to  have  the  beams  squally  distributed 
about  the  array  borealght;  N/2  beans  to  the  left  of  boreslght  and  N/2  beaji!:*  to 
the  right  of  boreslght.  (Boreslght  Is  defined  as  the  direction  perpendicular 
to  the  line  of  the  array.)  Also,  since  the  elements  are  equally  spaced,  a 
linear  phase  taper,  which  detersd-nas  a bean  direction.  Is  fomad  by  equal  phase 
Increments  between  elements  such  that  the  phase  of  the  k^*'  element  for  the  n^** 
beam  can  be  specified  as: 


^n,k  "^n,l 


where  4 . Is  determined  by  the  network  chosen.  Table  I shows  a Hating  of  the 

Ugl 

second  term  of  (1)  in  xultlplas  of  v for  an  elghc-elament  array.  The  synthanis 
of  the  network  required  to  form  those  phase  tapers  Is  described  elsewhere^  and 
Is  not  presented  here.  The  resultant  network  Is  shown  In  Figure  3. 

Figure  3 Indicates  a linear,  equally  spaced  array  of  antenna  elements. 

The  element  spacing,  however,  can  be  made  Irregular  to  sobm  extent^.  Figure  4 
shows  a plot  of  the  phase  taper  vs.  element  position  for  a typical  beam.  The 
plot  shows  the  phase  taper  for  the  eight  equally  spaced  elements  of  Figure  3a 
In  solid  line  and  the  taper  required  for  additional  elemants  If  the  array  were 
to  be  extended  without  disturbing  the  element  spacing  or  the  beam  direction  In 
dotted  Mne.  The  dashed  line  taper  would  be  required  If  the  elements  -4 
through  0 and  9 through  13  were  reversed  180*  In  phase.  This  can  generally  be 
accomplished  by  physically  rotating  the  antenna  elements  180*.  It  can  be 
seen  from  the  curve  tliat  since  the  phase  taper  can  be  drawn  as  a periodic 
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TABLE  I 

RELATIVE  PHASE  TERMS  - EIGHT  ELEMENT  LINEAR  ARRAY 


^SLtK.  (2ii-l)(k-l) 


IT 

M 

2 

3 

4 

5 

6 

7 

8 

1/8 

1/4 

3/8 

1/2 

5/8 

3/4 

7/8 

3/8 

3/4 

1-1/8 

1-1/2 

1-7/8 

2-1/4 

2-5/8 

5/8 

1-1/4 

1-7/8 

2-1/2 

3-1/8 

3-3/4 

4-3/8 

7/8 

1-3/4 

2-5/8 

3-1/2 

4-3/8 

5-1/4 

6-1/8 

1-1/8 

2-1/4 

3-3/8 

4-1/2 

5-5/8 

6-3/4 

7-7/8 

1-3/8 

2-3/4 

4-1/3 

5-1/2 

6-7/8 

8-1/4 

9-5/8 

1-5/8 

3-1/4 

4-7/8 

6-1/2 

8-1/8 

9-3/4 

11-3/8 

1-7/8 

3-3/4 

5-5/8 

7-1/2 

9-3/8 

11-1/4 

13-1/8 
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function,  elamants  within  the  bounds  of  the  original  array  can  be  moved  to 
positions  outside  these  bounds  where  the  required  phase,  including  the  180“ 
reversal,  is  no  different  from  that  at  the  original  element  position.  This  is 
done  with  elements  3 and  7 in  Figure  3b  to  produce  an  Irregularly  spaced  array. 

The  Butler  Beam  Forming  Network  can  then  be  used  for  both  equally  spaced 

linear  arrays  and  a class  of  randomly  spaced  linear  arrays.  It  is,  however, 

restricted  to  arrays  which  have  the  number  of  elements  equal  to  an  Integral 
2 

power  of  2.  Nolan  removed  this  restriction. 


IV.  THE  NOLAH  SYNTHESIS 


The  Btitler  Beam  Forming  Network  Is  reatrlcted  to  arrays  which  have  the 

number  of  elements  equal  to  an  Integral  power  of  2.  The  use  of  a Butler 

network  with  arrays  with  different  numbers  of  elements  usually  requires  the 

Introduction  of  some  loss  Into  the  network  or  the  use  of  couplers  with  more 

than  four  ports.  In  the  latter  case.  If  six-port  couplers  ware  utilised 

rather  than  four-port  couplers,  the  arrays  which  have  the  number  of  elements 

equal  to  an  Integral  power  of  three  could  be  utilized.  For  the  former  case,  an 

N element  array  could  be  used  with  an  M element  Butler  network  (M  greater  than 

or  equal  to  N)  with  the  unused  outputs  of  the  network  terminated.  This  results 

N 2 

in  a reduction  In  gain  by  a factor  (^)  . 


It  Is  possible,  however,  to  determine  a beam  forming  network  for  any 
arbitrary  number  of  array  elements  provided  that  the  network  transfer  function, 
A,  has  specific  properties.  In  general,  A can  be  written  us  an  NxN  matrix 
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*N1 


‘12 


‘22 


‘IN 


‘2N 


“N2 


• • 


*NN 


where  a , Is  the  transfer  coefficient  between  the  n beam  port  and  the  k 
n,k 

array  element  and 

a , ■ b , • where  b , Is  a real  number.  (3) 

n,k  n,k  n,k 


If  the  beam  forming  network  Is  to  be  lossless,  linear,  reciprocal,  and  matched, 
the  matrix  A must  have  the  following  properties: 


iO 


1. 


The  sum  of  the  squares  of  the  magnitudes  of  the  terms  In  any  row  or 
colunn  Is  unity. 

2.  The  sum  of  the  products  cf  the  terms  In  any  row  or  column  with  the 
conjugate  of  the  corresponding  terms  in  any  other  row  or  column  is 
zero. 

These  properties  are  an  expression  of  conservation  of  energy,  i.e.,  in  a 
lossless  matched  system,  the  power  in  equals  the  powez  out.  The  first  of  these 
simply  states  that  the  power  radiated  by  each  array  element  for  a given  beam  Is 
equal  to  the  power  input  to  the  corresponding  beam  port.  The  second  states 
that  If  a received  signal  Induces  voltages  In  the  array  elements  which  are 
equal  to  the  conjugate  of  the  transfer  coefficients  for  a given  beam,  then  all 
of  the  received  power  will  be  coupled  to  that  beam  port,  and  no  power  will 
couple  to  any  other  port.  Because  the  network  is  reciprocal,  It  may  be  reversed, 
the  beam  ports  and  element  ports  interchanged,  without  affecting  Its  perfor- 
mance; hence  the  statements  apply  to  both  rows  and  columns  of  the  transfer 
matrix. 

These  are  the  properties  of  a unitary  matrix  and  can  be  stated  mathema- 
tically In  matrix  notation  as 

a'^A  - I (4) 

where  I is  the  Identity  matrix  defined  by 

. ^ ^ n*k 

n,k"  1 Q n^k 

T 

and  ( ) represents  the  conjugate  transpose  of  the  matrix. 

One  of  the  properties  of  a unitary  matrix  Is  that  It  can  he  factored  Into 
the  product  of  elementary  unitary  matrices.  A particular  form  of  elementary 
unlt,*try  matrix  which  is  of  Interest  because  it  represents  a simple  r.f.  circuit 
Is 
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Fig.  5.  Circuit  representation  of  N x N elementary  unitary  matrix. 


Where  the  element  6 Is  a directional  coupler  with  the  properties  as  shown  In 
Figure  6. 


Fig.  6.  Di -actional  coupler  transfer  properties. 


This  elementary  circuit  provides  a means  for  determining  a beam  forming  network 
for  the  general  N element  array.  The  transfer  function  of  the  array  is  speci- 
fied by  (2)  and  (3)  and  can  be  drawn  In  circuit  form  as  shown  In  Figure  7. 


Fig.  7.  Beam  forming  network  before  reduction. 


The  circuit  A can  be  also  configured,  usiug  the  circuit  in  Figure  5,  ae 
shown  In  Figure  S: 


is-e-»a7Bi] 


Fig.  8.  Beam  forming  network  - first  Nolan  reduction. 

Comparing  Figures  7 and  8,  and  making  use  of  the  relations  In  Figure  6,  it  can 
be  seen  by  Inspection  that 

a , ■ a'  , cos  e e ••  a*  , sin  6 e (6) 

n,l  n,l  n»2 

-J-  j- 

a „ - a'  , sin  8 e ^ + a'  , cos  0 e ^ (7) 

or.  In  matrix  notation 
A'Ej^  - A 

The  circuit  A*  can  then  be  reduced  in  the  same  manner  as  was  the  circuit  A by 
oneratlng  on  terminals  1*  and  3 with  a similar  coupler,  such  that  the  circuit 
If.  Figure  9 is  obtained. 
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Continuing  the  same  process  for  N'l  couplers  results  in  the  circuit  In  Figure  10 


llt-t-uwl 


Fig.  10.  beam  forming  network  - N-lst  Nolan  reduction. 


N * 

where  the  transfer  function  a Is  an  an  N-1  matrix*  a matrix  of  order  less 

N* 

than  the  original  matrix  A.  The  circuit  A can  then  i.f:  reduced  In  the 
same  manner  as  was  A,  starting  with  the  2*  and  3*  ports  and  so  on*  until  the 
circuit  Is  completely  reduced  to  the  form  shown  In  Figure  11. 


) 2 3 4 • • • N -^element  PORTS 
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It  remains  now  to  detarmlne  the  values  of  6 and  ^ throughout  the  circuit, 
to  complete  the  synthesis.  The  two  preceding  circuits  show  that  there  Is  no 

i|i 

coupling  between  been  port  I and  the  port  2*  of  A end  simllerly  of  A*  because, 
as  shown  In  Figure  6,  an  input  from  the  right  of  any  coupler  Is  not  coupled  to 
its  bottom  port.  This  means  that  the  matrix  element  a*^  . la  equal  to  aero. 

This  allows  Equations  (6)  and  (7)  tc  be  solved  for  6^^  and  and  for  a^  j . 
Simllarlv,  a"j^  ^ Is  equal  to  zero  and  62  ^2  *^*^^™^**^*  l^ls 

procedure  can  be  continued,  and  the  circuit  synthesis  then  is  thus  one  requiring 
successive  solutions  to  equations  of  the  form  of  (6)  and  (7).  Using  matrix 
notation,  this  method  can  be  succinctly  expressed.  Since  A Is  unitary.  It  can 
be  expressed  as  the  product  of  the  factors 

^ “ ®N(N-l)/2  ®N(N-l)/2  -1  ®3®2®1 


where  Is  an  elementary  matrix  of  the  general  form 


E - 


COS  6 e 
0 


0 sin  6 e 


j(*  +7) 

sin  6 e 0 


COB  6 e 


0 

0 


T T 

Then  successive  post  multiplications  of  A by  E^  , E2  , etc. , will  eventually 
reduce  A to  an  Identity  matrix.  The  procedure  is  to  choose  the  6 , ^ pare** 

H u 

meters  In  E^  through  so  as  to  force  successive  terms  in  the  upper  row  of 

the  product  matrix  to  zero  until  ell  of  the  terms  in  the  row  ere  zero  except 
the  first  term.  Because  the  matrix  is  unitary,  the  first  term  in  the  row  will 
be  unity,  and  all  other  terms  In  the  first  column  will  be  zero.  This  reduces 
the  N X N transfer  matrix  tc  an  N-1  x N>1  transfer  matrix.  The  second  row  is 
reduced  in  a similar  way  by  proper  choice  of  the  next  N-2  values  of  E , etc., 

n 

until  all  values  of  6^^  and  ere  determined. 
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The  circuit  synthesised  by  this  sMthod  usee  fixed  pheee  shifts  and  quad- 
reture  couplers  which,  In  general,  are  not  3 dB  couplers.  These  can  be  made 
using  3 dB  couplers  using  the  circuit  In  Figure  12,  however. 


[r8-6-l8785| 


♦■—H  9 


Fig.  12.  Directional  coupler  circuit  using  3 tld  couplers. 


The  transfer  function  la  then: 


1 I cos  6 e 


eln  9 e 


2 I sin  6 e 


cos  9 e 


as  Is  required  for  the  coupler  characteristic  chosen  for  the  synthesis. 

This  synthesis  procedure  is  applicable  to  an  array  of  any  number  of 
elements  In  any  desired  configuration.  The  only  restriction  Is  that  the  desired 
transfer  awtrlx  be  unitary.  Figure  13  shows  the  circuit  for  a 5 element  array 
synthesized  by  this  procedure  where  the  terms  of  the  transfer  matrix  take  the 
form 

^ J(n-l)<k-l)^ 

n,k  ^ 

for  the  n^^  beam  port  and  the  k^^  array  element.  This  network  can  be  used  for 


a five  element  linear  array  with  equally  spaced  elements  and  uniform  amplitude 
aperture  distribution. 


18 


V. 


A NETWORK  SIMPLIFICATION 


The  form  of  the  beam  port/eleawnt  port  transfer  coaf flclenta  for  the 
network  of  Figure  13  la  generally  applicable  to  all  linear  equally  spaced 
arrays  and,  as  will  be  described  In  a later  section,  to  some  planar  arrays  as 
well.  For  uniform  Illumination,  the  desired  network  transfer  coefficients  will 
usually  take  the  form 


*n,k  /n 


1(n-l)(k-l)^ 


(8) 


This  can  be  simplified  somewhat  by  noting  that  the  N-(k-2)sc  terms  In  each  row 
have  the  values 

^ J(n-l)(N-(k-l)]^ 

*n,[N-(k-2)] 

^ -J(n-l)(k-l)^  j(n-l)2it 

• — e e 

j<n-l)2it 

and  since  e >1 


*n,[n-(k-2)] 


^ -J(n-l)(k-l)-^j 

I — e 

/n 


which  Is  the  conjugate  of  the  term  a . . A similar  argument  for  the 

n,K 

N-(n-2)8t  terms  In  each  column  can  also  be  made.  This  suggests  a possible 
simplification  of  the  network  to  be  synthesised. 

Consider  a 3 dB^coupler  circuit  with  inputs  /2  cos  a and  /2  sin  a.  As 
Indicated  In  Figure  14,  the  outputs  will  take  the  form  e^'^  and  e 
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Fig.  13.  Five-element  array  BFN 


Figure  IS  shows  that  the  matrix  B can  be  treated  as  two  separate  matrices 

N+1  ^ 

and  B2  where  B^  Is  comprised  of  the  first  columns  of  B,  and  B^  is  comprised 
of  the  remaining  columns. 
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Examination  of  B^  shows  that  the  n row  Is  ids  ' to  the  N-(n-2)st  row 

since  cos  a « cosC*^).  Examination  of  B,  shows  thk...  the  N-(n-2)st  row  Is  the 
th  ^ 

negative  of  the  n row  since  sin  a > -sin(-a).  This  suggests  that  the  number 
of  rows  In  each  matrix  might  be  reduced  since  the  difference  between  the  for- 
matlon  of  the  n beam  and  the  (N-(n-2)st  beam  Is  simply  the  sign  of  the  terms 
in  the  B^  matrix. 
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A coupler  with  the  properties  shown  In  Figure  16  selves  the  purpose  for  this 


Fig.  16.  3 dB  coupler  circuit  transfer  properties. 


and  the  network  becomes  that  shown  In  Figure  17.  The  90**  phase  shifter  at 
the  beam  Inputs  Indicated  In  Figure  16  Is  omitted  In  Figure  17  since  this 
affects  only  the  relative  phase  of  the  beam.  The  transfer  matrices  of  the 


networks  and  C2  are 
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In  general,  the  Nolan  synthesis  of  an  N x N network  will  require  — — 
of  the  6 type  couplers  described  In  the  previous  section  or,  as  previously 
Indicated,  twice  as  many  3 dB  couplers,  NCN-l).  The  five  element  array  in 
Figure  13  uses  10  of  the  6 couplers,  and  If  fabricated  with  3 dB  couplers  would 
require  20  of  these.  Examination  of  Figure  13  shows  that  this  number  can  be 
reduced  to  18  since  the  two  45°  couplers  are  already  3 dB  couplers.  The  pre- 
ceding simplification  for  transfer  coefficients  of  the  form  (8)  where  N Is  odd 

j .V  ..u  4 r . ..  u j ..V  N-1  N-1 

require  tV»e  synthesis  of  two  networks,  one  x -y-  and  the  other  — y x ■~y» 

Together  these  matrices  would  require  of  the  6 type  couplers  and  N-1 

N^-1  ^ X / 

3 dB  couplers  or  - — 3 dB  couplers.  The  total  number  of  3 dB  couplers  required 
is  reduced  by  a factor  "yy*  This  reduces  the  number  required  for  a five- 
element  array  to  twelve  3 dB  couplers.  The  five-element  array  network  has  been 
synthesized  using  the  simplification,  and  the  result  Is  shown  In  Figure  18. 
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Thus  far,  the  discussion  of  simplification  of  the  network  has  been  con- 
cerned with  networks  for  arrays  which  consist  of  an  cdd  number  of  antenna  ele- 
me-  Because  of  this,  the  first  column  of  the  A transfer  matrix  was  left 
unu  arbed  In  the  conversion  to  the  B transfer  matrix.  If  the  number  of 
elem  s and  beams  Is  even,  a matrix  with  transfer  coefficients  of  the  form  (8) 
is  generally  unsuitable.  This  matrix  would  Include  an  eleme.nt-to-element  phase 
shift  xBO”.  This  results  In  a beam,  for  a linear  array,  which  points  along 
the  axis  of  the  array  In  both  directions.  A more  suitable  transfer  matrix  for 
even  n.  ./ered  arrays  Is  the  form  used  for  the  Butler  Matrix 

(9) 


^ j(2n-l)(k-l)-^ 

a , •=  — e 
n,k 


/n 

This  form  results  from  adding  a phase  shift  of  (k-1)^  to  the  k^^  element  port 
of  a matrix  with  transfer  coefficients  of  the  form  (8). 

Using  a desired  transfer  matrix  of  the  form  (9)  with  N even,  the  rows 
and  columns  of  the  matrix  having  a 0-180°  or  + 90°  relationship  may  first  be 


This  reduces  the  transfer  matrix  to  two  ^ ^ Y 


combined  with  3 dB  couplers, 
transfer  matrices.  If  ^ odd,  then  each  matrix  may  be  simplified  as  descriued 
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! 1 


i I 


for  arrays  with  odd  numbers  of  elements.  If  ^ Is  even,  then  the  matrix  may  be 
further  reduced.  The  simplification  of  the  eight  element  network  la  presented 
as  an  example. 

Using  (9) , the  A matrix  Is 
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Columns  l&S,  2&6,  3&7,  and  4 & 8 have  ^ ± ^ phase  relationship  on  a term  by 
term  basis.  These  are  combined  using  3 dB  couplers  to  form  the  circuit  In 
Figure  19. 
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The  end  B2  transfer  natrices  becone: 
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Note  now  that  adding  ir/A  to  columns  2 and  4 of  B^  and  to  columns  1 and  3 of 
again  produce  a >(if/2)  relationship  as  before.  Therefore  using  a coupler 
which  results  in  a +(if/2)  relationship  and  subtracting  w/4  produces  .he 
desired  result.  Figure  20  shows  the  resultant  circuit. 
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The  C matrices  are 
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And  note  here  that  the  addition  of  3ff/8  to  column  2 of  and  column  1 of 
results  In  a quadrature  relationship.  Also  the  addition  of  ir/8  to 
column  1 of  and  column  2 of  producea  a quadrature  relationship. 
Subractlng  theiie,  as  before,  from  the  coupler  outputs  producss  the  circuit, 
neglecting  tcniis  which  only  affect  the  relative  phase  ef  the  bens,  shown 
in  Figure  21.  TTils  circuit  Is  Identical  to  that  shown  it  Figure  3 for  the 
eight  element  Butler  beam- forming  network. 

f 
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It  should  be  noted  that  1£>  during  the  simplifying  process,  the  required 
transfer  matrices  had  an  odd  number  of  rows  and  columns,  as  would  be  the 
case  for  a 6 element  network,  the  simplification  procedure  described  for 
odd  numbered  arrays  Is  applicable. 

The  result  In  the  preceding  example  Indicates  that  the  Butler  matrix  Is 
a special  case  of  a more  general  beam  forming  network  synthesis  procedure. 

It  Is  possible  to  apply  the  Nolan  synthesis  p. jcedure  directly  without  using 
the  procedures  described  In  this  section  and  obtain  the  same  result.  This 
P ocedure,  however.  Involves  reordering  the  columns  of  the  transfer  matrix 
after  each  operation.  The  methods  presented  here  for  simplification  of  the 
networks  involve  less  effort  than  does  the  direct  application  of  the  Nolan 
procedure. 
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VI.  PLANAR  ARRAYS 


The  phase  taper  for  a linear  array  can  be  represented  graphically  as  a 
line  as  In  Figure  4.  The  phase  taper  for  a planar  array  tikes  on  an  added 
dimension  and  can  be  represented  as  a plane.  If  the  phase  argument  for  an 
N element  linear  array  Is  as  In  (5)  (l.e.,  (n-1) (k-1)-^) , then  the  phase  taper 
will  begin  a new  cycle  at  the  N-t-lst  element.  Application  of  this  property  to 
planar  arrays  results  In  the  synthesis  of  multlbeam  planar  arrays  described 
by  Shelton. First,  the  geometric  configuration  of  the  desired  planar  array 
Is  selected  and  a reference  element  (k-1)  Is  arbitrarily  chosen  as  In  the  seven 
element  hexagon  In  Figure  22. 
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REFi'RErCE *-0  O O 
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Fig.  22.  Seven-element  array. 


Next,  this  array  Is  Inserted  Into  a lattice  structure  of  "phantom"  arrays  as  In 
Figure  23. 
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It  can  then  be  seen  that  the  phase  of  the  reference  element  of  arrays  A&B 
must  differ  from  the  reference  element  of  the  "real"  array  by  In/n  for  the 
beam.  Since  the  phase  taper  Is  planar,  one  can  write 

m X + m Y - 2tt  (n-1)  (10) 

X y 

where  m and  are  the  slopes  of  the  phase  taper  In  the  X and  Y directions, 

^ I 

respectively.  If  X and  Y are  measured  In  terms  of  the  number  of  spaces  between 
elements,  then  Equation  (10)  can  be  written  as 

K^(n-l)  ^ X + Ky(n-l)  Y - 2n  (r  (11) 


where  and  are  Integers.  The  equation  simplifies  to 

X X + K Y • N (12) 

X y 

and  for  the  example  of  the  7-alewent  array  with  the  "phantom"  arrays  A and  B, 
the  two  equations 


3K  - K - 7 

X y 


K + 2K  - 7 

X y 


can  be  written  for  which 


K - 3 

X 


K - 2. 

y 

This  determines  the  phase  Increments  for  the  7-element  array.  These  are 
tabulated  for  k»l  through  7 In  Figure  24. 
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Fig.  24.  Flement  Phases  as  a multiple  of  ^ for  a 7-element  hexagonal 
array. 
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Note  that  If  these  phase  Increments  are  tabulated  with  element  numbers  assigned 
as  in  Figure  25,  the  desired  transfer  matrix  becomes: 


A 


(1//7)  e^<2r/7) 


0 

0 

0 

0 

0 


t 0 
I 

Li> 


0 

1 

2 

3 

4 

5 

6 


0 

2 

4 
6 
1 
3 

5 


0 

3 
6 
2 
5 
1 

4 


0 0 0 n-1 

4 5 6 n-2 

1 3 5 n-3 

5 1 4 n-4 

2 6 3 n-5 

6 4 2 n-6 

3 2 1 n-7j 


3 e 

O O 


'O  «0  TO 


20  so 

Fig.  25.  Seven-element  hexagonal  array. 


This  matrix  has  transi ar  coefficients  of  the  form 


a 


n,k 


j(n-l)(k-l)-^ 
(1//7)  e ” 


which  is  no  different  from  the  matrix  previously  described  for  a linear  array. 
The  beam  forming  network  for  this  type  of  array  is  identical  to  that  used  for 
a linear  array. 
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Although  this  procedure  allows  the  beam  forming  network  to  be  applied 
to  planar  arrays,  It  Is  somewhat  restricted.  First,  the  array  elements  must  be 
arranged  such  that  they  are  on  a regular  lattice  such  as  in  the  example; 
second,  not  all  arrays  even  it  configured  in  this  way  will  result  in  use  of  the 
N beam-forming  network.  The  8-element  array,  for  example,  configured  in  a 
rectangular  lattice  reduces  to  two  four  element  arrays  by  this  method. 

Random  element  spacing,  subject  to  the  restriction  that  the  array  elements 
lit!  on  the  lattice,  is  possible.  This  Is  an  extension  of  the  principle  described 
for  the  Butler  array.  The  appropriate  elements  in  the  "real"  array  are  replaced 
by  the  elements  in  the  "phantom"  array. 
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VII.  THE  GENERAL  ARRAY  CASE 


In  the  beginning  of  the  discussion,  the  restrictions  placed  on  the  antenna 
array  used  with  the  beam  forming  network,  were  that  It  was  required  to  be  a 
linear,  equally  spaced  array  of  similar  antenna  elements  of  number  equal  to  an 
Integral  power  of  2.  Subsequent  discussion  removed  the  equal  spacing  require- 
ment and  substituted  the  requirement  that  the  element  spaclngs  must  be  Integral 
multiples  of  some  arbltarlly  chosen  basic  spacing  Interval.  Section  IV  extends 
the  array  to  allow  any  number  of  antenna  elements  by  using  the  Nolw  synthesis 
procedure.  The  restriction  that  the  array  be  a linear  array  was  removed  In 
Section  VI  and  the  use  of  planar  arrays  configured  on  a regular  lattice  was 
discussed. 

There  are.  Inevitably,  some  design  situations  where  more  general  types  of 
arrays  are  required.  For  example,  it  may  be  necessary  to  configure  the  array 
In  three  dimensions  rather  chan  one  or  two;  or,  the  array  element  patterns  may 
be  slfjnlf icantly  different  for  one  reason  or  another. 

T'-  N-tlan  oypthesls  Is  confined  by  only  one  restriction,  that  the  speci- 
fied tri'  .rer  fimctlon  be  a unitary  matrix.  This  is  another  way  of 

saylt..,  hat  the  beam-forming  network  is  lossless.  Throughout  the  discussion 
thus  far,  only  two  unitary  transfer  matrices  have  been  discussed,  those  of 
Equations  (8)  and  (9).  These  do  not  exhaust  the  possibilities,  however. 

In  the  general  case,  a transfer  matrix  which  Is  not  unitary  may  be  desired 
In  order  to  achieve  a particular  set  of  array  excitations.  Such  a network 
will  be  lossy.  It  is  d-*''  'ult  to  determine  a minimum  loss  network  which  will 
have  the  des.ii  : cransf  . junction;  and  If  the  network  can  be  determined.  It  Is 
very  likely  to  be  too  lossy  and  to  produce  beams  of  unequal  gain.  Usually, 
however,  there  Is  some  tolerance  on  the  desired  beam  directions  so  that  the 
desired  transfer  function  )e  considered  to  be  approximate.  Allowing  some 
deviation  from  the  desire  ..ansfer  function,  the  problem  then  becomes  one  of 

finding  a unitary  transfer  function  which  best  approximates  In  some  sense  the 
desired  transfer  function. 
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One  methoJ  of  appro:  imatlon  Is  to  minimize  the  mean  square  difference 
between  the  terms  of  the  unitary  matrix  and  the  deelred  network  transfer 
matrix.  Using  this  as  a basis  for  approximation,  the  unitary  matrix  which  best 
approximates  a given  desired  transfer  function  Is  given  as 
T 

A - VT 

where  A Is  the  resultant  unitary  matrix  and  V and  T ara  eigenvector  matrices 

T T 

associated  with  the  product  matrices  SS  and  S respectively . S Is  the 

T 

desired  transfer  matrix  and  ( ) indicates  the  conjugate  transpose  of  a matrix. 
The  derivation  of  this  expression  Is  described  In  Appendix  I.  The  matrix  A 
will  not  necessarily  represent  a transfer  function  which  will  produce  accep- 
table antenna  bears.  Much  depends  on  the  original  desired  transfer  function  S. 
If  the  desired  beams  are  uniformly  distributed  in  the  region  between  grating 
lobe  positions  of  the  array,  t:  an  A will  provide  a good  approximation  to  the 
desired  signal  matrix.  If,  however,  the  desired  beam^-  are  closely  grouped  with 
a relatively  large  angular  separation  between  this  group  and  the  grating  lobe 
positions,  then  the  approximation  will  probably  be  a poor  one.  In  more  rigorous 
terms,  the  xoore  completely  the  desired  beam  vectors  span  tho  N dimensional 
vector  space  of  the  NxN  desired  transfer  matrix,  the  closer  the  approximation 
of  A to  S.  This  Is  simply  a result  of  the  fact  that  beam  vectors  of  A will 
always  completely  span  the  vector  space. 

The  antecna  beams  obtained  with  a beam-forming  network  A will  usually  be 
somewhat  different  in  gain  with  respect  to  one  another,  the  differences  depend- 
ing on  the  closeness  of  the  approximation  to  the  transfer  function  S.  (This 
assumes  that  S specifies  equal  gain  beams  although  this  is  not  necessary.  If  S 
specifies  unequal  gains  than  the  statement  applies  In  a relative  sense.)  This 
results  from  the  fact  that  the  sidelobe  structures  of  the  individual  beams  jflll 
differ  and  that  the  beamwldths  and  shapes  of  the  beams  will  also  not  be  Iden- 
tical. 

The  approximation  described  coupled  with  the  Nolan  synthesis  pro/idos  some 
significant  advantages  in  most  cases,  however.  Together,  they  represent  a 
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straightforward  procedure  for  the  design  of  losslees  beam-fomlng  networks  for 
arbitrary  arrays.  Since  the  desired  transfer  matrix  S is  completely  arbitrary 
the  array  can  take  any  conceivable  configuration.  The  array  element  spacing 
can  be  anything.  The  array  can  be  linear,  planar,  or  three  dimensional;  it  can 
consist  of  regularly  spaced  or  randomly  spaced  elements.  The  array  can  be  a 
filled  array  or  a thinned  array,  and  it  can  even  be  "tade  up  of  elements  with 
differing  pattern  gain  characteristics.  In  addition,  it  is  possible  to  specify 
aperture  distributions  to  produce  low  sidelobes  or  to  shape  the  beam.  The 
degree  to  which  the  resultant  design  satisfies  these  requirements  is,  in 
general,  dependent  only  upon  how  completely  the  specified  signal  vectors  span 
the  N dimensional  vector  space  of  the  matrix. 

The  results  of  this  approximation  have  bean  applied  to  the  synthesis  of 
beam-forming  networks  for  randomly  spaced  thinned  planar  arrays  of  similar 
elements.  The  characteristics  of  two  such  arrays  were  computed.  The  first,  a 7- 
element  array  with  a single  central  element  and  6 elements  randomly  spaced  on 
the  circumference  of  10  foot  radius  circle.  The  desired  beams  were  specified 
to  be  of  equal  gain  and  to  be  located  on  the  vertices  of  a hexagon  at  an 
angular  separation  of  7°  from  the  array  boreslght.  A single  beam  was  specified 
to  be  on  boresight.  Figures  26  and  27  show  the  computed  1 dB  and  3 dB  contours 
for  the  beams  formed  at  350  MHz  using  the  transfer  function  A. 

The  second  array  computed  was  a 19-element  array  with  3 elements  randomly 
located  on  a 15  foot  diameter  circle,  6 elements  on  30  foot  diameter  circle, 
and  10  elements  on  a 45  foot  diameter  circle.  The  de  red  beam  directions  were 
again  chosen  in  a hexagonal  partem  to  form  two  concent  c hexagons  with  a beam 
on  the  rray  boresight.  Figures  28  and  29  show  the  computed  1 dB  contours  at 
350  MHz  for  this  array  with  a beam-forming  network  with  the  unitary  transfer 
function  A and  the  array  configuration,  respectively. 
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Fig.  26.  Random  circular  array  1 dB  contours. 
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Fig.  26.  19>b«m  array. 
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VIII.  FUMMARY 


This  discussion  has  attempted  to  describe  possible  methods  of  deriving 
lossless  beam-forming  networks  for  various  types  of  arrays  beginning  with  the 
well-known  Butler  array  and  extending  this  array  to  unequally  spaced  linear 
arrays,  planar  arrays  and  finally  to  the  completely  general  type  of  array.  For 
this  last,  it  is  indicated  that  it  is  not  always  possible  to  synthesize  a 
lossless  BFN  but  that  there  exists  a least  mean  square  optimization  which  leads 
to  the  synthesis  of  a lossless  BFN  which  provides  perfo  mar.ce  which  approxi- 
mates the  desired  performance. 
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APPENDIX  I 


I 


Denoting  the  optimum  unitary  transfer  function  by  A then  as  In  Equation  (4) 
T 

A A - I 

T 

and  denoting  the  desired  transfer  function  by  S and  premultiplying  by  A 
results  lu 


A^S  - B (A-1) 

where  B Is  approximately  equal  to  1.  If  S Is  represented  as  the  sum  of  the 
matrix  A and  some  matrix  R such  that 

A + R » S (A-2) 


then 


A^S  - A^(A  + F)  - I + A^R  - B 


and 


A^R  - B - 1 

which  Indicates  that  If  A Is  a non  trlval  (not  all  zeros)  matrix,  that  It  is 
desirable  to  minimize  R In  order  to  make  B approximate  I as  closely  as  pos- 
sible. 

One  way  of  minimizing  R Is  to  minimize  the  sum  of  the  squares  of  the 
magnitudes  of  each  term  In  R.  Or, 
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i 

i 


or,  with  more  convenient  notation 

I I r|  I ^ ■ minimum  (A-3) 

In  order  to  determine  the  optimum  transfer  matrix  A which  will  minimize  R, 

rgi 

It  Is  necessary  to  assert  two  matrix  transformations.  ''  First,  the  eigenvalue 
transformat Ion 

- A (A-4) 

where  M Is  any  matrix  with  n rows  and  P columns  and  la  an  NxN  unitary 

matrix.  The  matrix  A Is  diagonal  with  each  term  on  the  diagonal  a real  number 

T 

called  an  eigenvalue  of  the  matrix  product  MM  . The  columns  of  are  called 

T T 

the  eigenvectors  of  the  matrix  MM  . If  the  matrix  MM  Is  known,  then  the 
eigenvalues  A^  and  the  eigenvectors  of  can  then  be  determined  by  straight- 
forward (but  tedious)  computational  methods.  Second,  using  the  notation  of 


(A-3) 

I ( U2NI  I ^ - I I N U3I  I ^ - I I N|  1 ^ (A-5) 

where  and  arc  any  unitary  NxN  matrices  and  N Is  any  NxN  matrix.  This  can 
be  demonstrated  by  carrying  out  the  summation  on  a term  by  term  basis  and  using 
the  properties  of  the  unitary  matrix. 


2 

Referring  back  to  Equation  (A-3) , it  Is  desired  to  minimize  1 1 R|  | and 
hence  from  Equation  (A-2) 


[ ! 


There  is  a matrix  transformation  which  v?lll  diagonalize  the  N by  N matrix  S 
such  that 


V^ST  - I (A.7) 

Where  V and  T are  unitary  matricea  and  Z is  a diagonal  matrix  with  positive 
real  terms  on  the  diagonal.  This  is  shown  by  premult ip lying  (A-7)  by  V and 
postmultlplylng  by  and 

T T T 

W STT  - VZT 

and  since 

T T 

W - TT^  - I 


VZT 


From  (A-8) , S ■ TZV^,  since  Z is  diagonal  positive  real  and 


then 


SS^  - V Z t\ 


2 T 
« V E^ 


and 


T T 2 

V SS  V - E 
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“7 


Since  the  matrix  Is  diagonal  positive  real,  then  from  Equation  (A~4),  V Is 

T 

the  eigenvector  matrix  of  SS  and  the  terms  of  the  matrix  Z,  o^,  are  the  square 
roots  of  the  eigenvalues  Similarly, 


T 

S S 


T T 

T Z V V Z T 


2 T 

- T Z T 


and 


T T 

T S ST 


where  it  is  seen  that  T is  the  eigenvector  matrix  of  S S. 
Using  Eq.  (A5)  to  re-write  Equation  (A-6)  provides 


s-a| 


V^ST  - v”at| 


mlr.lm>ia 


or 


z - v'^at!  I ^ 


minimum 


The  probleo'.  then  is  to  minimize  | | Z-P  | ^ where  Z Is  diagonal  positive  real  and 
P is  the  matrix  V^AT. 

Writing  as  a term  by  term  summation 

, N K - 

2 _ ^ ^ I ^ „l  2 


r-p 


N K 

Z Z I Z-P[ 
1-1  1-1  ij 


but  since  Z Is  diagonal 


E-P|| 


H N 


N 


z z ip.,r+  z io^-p^,i 


1-1  j-1 
i'/j 


ij' 


1-1 
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and  expanding 


N N 

" E 
i-.  J-1 

ir‘J 


N 

+ Z 


i-1 


+ 


+ 0 


where  ( ) indicates  the  conjugate. 
And  rewriting 


E-P 


N N 
E E 
1-1  j-1 


'IJ' 


N 

+ E 
1-1 


” 2 

Now  the  double  summation  for  the  first  term  Is  positive  and  also  E a.  Is 

1-1 

some  positive  number  which  la  dependent  upon  how  the  original  desired  transfer 

2 

matrix  was  determined.  Sin '.e  1 1 E-P|  | must  be  greater  than  or  equal  to  zero, 
It  Is  minimized  when  the  negative  term  of  Equation  (A-9)  Is  maximized,  l.e.. 


" * 

^ <Ji(Pil  P (A-10) 

1-1 

Nov  the  matrix  P is  unitary  since 

P P - T^A^W  AT  - T A AT  - T T - I 

and  the  maximum  possible  value  of  Equation  (A-10)  Is  achieved  when  p^^-l'^jO. 
If  this  Is  true,  then  the  matrix  P is  the  identity  matrix  I and 


V^AT  - I 


from  which  Is  follows  that  the  optimum  unitary  beam  forming  network  transfer 
function  Is 
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where  V and  T are  the  elganvactor  wtrlcaa  aaaociatad  with  tha  product  matrices 
SS  and  S S of  the  desired  non-ualtary  transfer  function  S. 


I 
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